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ABSTRACT 

Fermi acceleration can develop efficiently at relativistic collisionless shock waves provided 
the upstream (unshocked) plasma is weakly magnetized. This has been both indicated by an- 
alytical theory and observed in numerical particle-in-cell simulations. At low magnetization, 
the large size of the shock precursor indeed provides enough time for electromagnetic micro- 
instabilities to grow and such micro-instabilities generate small scale turbulence that in turn 
provides the scattering required for particles to undergo Fermi cycles at superluminal rela- 
tivistic shock waves. The present paper extends our previous analysis on the development of 
these micro-instabilities to account for the finite angular dispersion of the beam of reflected 
and accelerated particles and to account for the expected heating of the upstream electrons 
in the shock precursor Indeed we argue that the electrons can be significantly heated during 
their travel in the shock precursor and that they may even reach equipartition with protons, in 
agreement with recent numerical simulations. We show that the oblique two stream instability 
may operate down to values of the shock Lorentz factor ^ 10 (corresponding to a rela- 
tively large angular dispersion of the beam) as long as the electrons of the upstream plasma 
remain cold, while the filamentation instability is strongly inhibited in this limit; however, as 
electrons get heated to relativistic temperatures, the situation becomes opposite and the two 
stream instability becomes inhibited while the filamentation mode becomes efficient, even 
at moderate values of the shock Lorentz factor The peak wavelength of these instabilities 
migrates from the inertial electron scale towards the proton inertial scale as the background 
electrons get progressively heated during the crossing of the shock precursor We also discuss 
the emergence and the role of current driven instabilities upstream of the shock. In particular, 
we show that the returning and accelerated particles give rise to a transverse current through 
their rotation in the background magnetic field. We find that the compensating current in the 
background plasma can lead to a Buneman instability which provides an efficient source of 
electron heating. 
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1 INTRODUCTION 

The acceleration of particles through repeated interactions with the 
electromagnetic fields up- and down-stream of a collisionless shock 
front is generally taken as the source the non-thermal particle pop- 
ulations that are detected in powerful astrophysical outflows. The 
physics of this Fermi mechanism is relatively well understood for 
non-relativistic shock velocities, at least in the test particle limit in 
which one neglects the backreaction of the accelerated particles on 
the shock environment. In the relativistic limit however, in which 
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the shock moves toward the upstream (unshocked) medium with 
a bulk Lorentz factor 2> 1, the situation becomes more in- 
tricate, mostly because the shock now moves about as fast as the 
accelerated particle. In non-relativistic shocks one can distinguish 
the details of the shock structure from that of the Fermi process, as 
this latter takes place on spatial scales much larger than the shock 
thickness. On the contrary, in relativistic shocks, the structure of 
the shock plays a significant role in the Fermi process and in the 
process of generation of the electromagnetic turbulence. One im- 
portant issue in particular, is the reflection of a part of the cold up- 
stream particles at the shock front, which together with the Fermi 
accelerated particles, travel back upstream and initiate instabilities 
in the upstream plasma. 
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These instabilities play a key role in the development of the 
Fermi process and consequently, in the radiative signatures of rel- 
ativistic outflows. Indeed, at magnetized oblique ultra-relativistic 
shock waves, the Fermi process may develop only if sufficient 
micro-turbulence has been generated in order to prevent the advec- 
tion of the accelerated particles away from the shock front into the 
far downstream plasma (Begelman & Kirk 1990, Niemiec & Os- 
trowski 2006, Lemoine et al. 2006). And, for the vast majority of 
them, ultra-relativistic shock waves are oblique: the opposite paral- 
lel configuration only takes effect when the magnetic field lines are 
oriented to an angle with the shock normal smaller than 1 / 7sh (in 
the unshocked plasma or upstream rest frame). As to the meaning 
of "magnetized" shock waves, it is slightly ambiguous and context 
dependent; above, it is meant: with a pre-existing (upstream) mag- 
netic field Bu such that, after compression through the shock and 
Lorentz transform into the downstream frame, the typical Larmor 
radius of accelerated particles is much smaller than the width of the 
blast, which is of the order of R/'y^ti in this frame (i? denotes the 
shock radius). For smaller values of B^, the accelerated particles 
can explore the whole blast and beyond, hence they may scatter on 
magnetic inhomogeneities sourced elsewhere, which would allow 
their repeated interactions with the shock. Nevertheless, for such 
low values of B^, the generation of micro-turbulence is guaran- 
teed, as will be discussed in the following, hence the Fermi pro- 
cess should develop independently of any extra source of magnetic 
inhomogeneitie|3 Furthermore, there exist evidence for amplifica- 
tion of the magnetic field upstream of a relativistic shock, see e.g. 
Li & Waxman (2006), Li & Zhao (2011), which can be seen as a 
hallmark of the development of micro-instabilities upstream of the 
shock. 

On microscopic plasma scales, several instabilities may de- 
velop in the precursor of the shock front through the penetration of 
the beam of returning and accelerated particles into the upstream 
plasma. The filamentation instability in particular - often referred 
to as the Weibel instability in the literature - has been proposed 
as the agent of amplification of the magnetic field to the level in- 
ferred from the synchrotron interpretation of gamma-ray burst af- 
terglows (Gruzinov & Waxman 1999, Medvedev & Loeb 1999; 
see also Wiersma & Achterberg 2004, Lyubarsky & Eichler 2006, 
Achterberg & Wiersma 2007, Achterberg et al. 2007). Recently, 
it has been found that the growth rate of two-stream instabilities 
in an oblique configuration (Bret, Firpo & Deutsch 2005) is larger 
than that of the filamentation mode (Bret 2009, Lemoine & Pel- 
letier 2010). Which instability grows faster is a key question here, 
as the relativistic velocity of the shock wave and the level of up- 
stream magnetization strongly limit the penetration length scale of 
the accelerated particles into the upstream plasma before these par- 
ticles are caught back by the shock wave (Milosavljevic & Nakar 
2006, Pelletier et al. 2009); and as mentioned before, if such in- 
stabilities cannot grow, then Fermi acceleration cannot develop (at 
least, in the absence of extra sources of turbulence). 

The detailed conditions under which the above instabilities 
and other relevant modes can grow, together with the develop- 

^ More generally speaking, Fermi acceleration requires micro-turbulence 
to be generated at distances smaller than r'L.o, with r'L.o the Larmor ra- 
dius in the downstream frame shock compressed background magnetic field 
(Pelletier et al. 2009). While this is guaranteed for micro-instabilities when 
they can be triggered in the precursor, this constraint can be turned into a 
lower bound on the growth rate of extra sources of turbulence, which of 
course differs from the constraint on the growth of micro-instabilities dis- 
cussed here and in Lemoine & Pelletier (2010). 



ment of Fermi cycles as a function of the shock velocity and the 
upstream magnetization have been discussed in a previous paper 
(Lemoine & Pelletier 2010). That paper assumed a charge neutral- 
ized beam without angular dispersion propagating in a cold back- 
ground plasma. It is the aim of the present paper to revisit these as- 
sumptions and to extend the previous calculations to a more generic 
situation, in which the beam angular dispersion (which is fixed by 
kinematics, as discussed below) is taken into account, and in which 
the effects of heating of the upstream plasma are considered. The 
emergence of current instabilities is also discussed. The present 
study is warranted in particular by recent work (Lyubarsky & Eich- 
ler 2006, Rabinak et al. 2011) that showed that the small opening 
angle of the beam suffices to prevent the onset of the filamenta- 
tion instability in electron-proton plasmas of shock Lorentz factor 
7sh ^ 100, independently of the length scale of the precursor. The 
present study recovers this result for a cold background plasma, 
but it also argues that: (i) the two-stream instability and more pre- 
cisely its oblique mode version is less sensitive to the beam angular 
dispersion as it can develop down to 7ah ~ 10 in the same condi- 
tions; (ii) the heating of the background plasma has a major effect 
in that it may help sustain the development of instabilities down to 
7sh ~ 10. 

The present paper is laid out as follows. In Section|2] we first 
review the current understanding on the structure of the shock, re- 
lying on the results of the most recent particle-in-cell (PIC) sim- 
ulations of relativistic collisionless shocks. We also discuss in a 
general way the instabilities that may develop as a function of the 
two main characteristics of the shock: upstream magnetization and 
shock velocity. In Section [3] we then discuss how the filamenta- 
tion and two-stream instabilities are affected by the finite angular 
dispersion of the beam of returning particles, assuming the back- 
ground plasma to remain cold. In Section|4]we generalize the calcu- 
lations of Sec.|3]to a background plasma composed of cold protons 
and relativistically hot electrons. Finally, in Sec.|5]we discuss the 
possibility of current instabilities and their role in shaping the pre- 
cursor. We summarize our results and draw conclusions in Sec. [6] 
In Appendix lAl we explicit the susceptibility tensor of the beam of 
finite angular dispersion, modeled with a waterbag distribution in 
the transverse momentum direction (transverse with respect to the 
shock normal) and a Dirac distribution in the parallel momentum 
direction. 



2 GENERAL CONSIDERATIONS 
2.1 Shock structure and precursor 

Let us first describe the general structure of a relativistic colli- 
sionless shock and introduce the main quantities of interest for the 
present discussion. 

The shock is characterized by a small number of parame- 
ters: (i) the composition ahead of the shock front, in the upstream 
plasma - in what follows, we assume except otherwise noted that 
the plasma is composed of electrons and protons with densities 
Ue — Up = Uu', (ii) the equation of state of the upstream plasma, 
which will be taken cold or composed of relativistically hot elec- 
trons but cold protons (see also further below for more details on 
this point); (iii) the shock velocity, written as /3shC as measured in 
the upstream frame, with corresponding Lorentz factor ^^h', (iv) the 
magnetization of the upstream plasma, written CTu and defined as 
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The magnetic field strength is here measured in the upstream 
frame. The shock structure also depends on the obliquity of the 
background magnetic field, but in what follows and except oth- 
erwise noted, we will focus on the generic superluminal case in 
which the angle with respect to the shock normal exceeds l/7sh. 
Then, in the shock front rest frame, the magnetic field lies mostly 
perpendicular to the flow direction, as the transverse components 
are amplified by a factor 7sh while the parallel component re- 
mains unchanged by the Lorentz transform. For an oblique shock 
wave, the above magnetization parameter thus measures (up to a 
sin^ 0B ~ 1 factor, with 0s the angle between the shock normal 
and the magnetic field in the upstream rest frame) the ratio of the 
magnetic energy in the shock front frame relative to the incoming 
matter flux ^%in^mp(? crossing the shock front. 

CoUisionless shock waves can be mediated by three generic 
types of mechanisms: by an electrostatic potential for incoming 
protons in an e — p shock, by a magnetic barrier or, if the mag- 
netic field is negligible, by a reflection due to the ponderomotive 
force of growing waves ahead of the shock transition layer. In the 
case of the external shock of a gamma-ray burst outflow propagat- 
ing in the interstellar mediiun, for example, the reflection can take 
place on a potential barrier or through a ponderomotive force, as the 
ambient magnetic field is then particularly weak, Uu ~ 10^^ for 
Bu ~ 1 and riu ~ 1 cm~^. However, even weak, the ambient 
magnetic field can play an important role, as will be seen further 
on. 

The magnetic field at a relativistic shock front is associated 
with a motional electric field E^^^Yi ~ Ah-Bu|sh when measured 
in the shock front frame (as indicated by the |sh subscript), with 
Su|gh — 7shSu (again, up to a factor sin 0b). Since the magnetic 
field is frozen in most parts of the plasma, the transverse component 
is further increased by the velocity decrease as the upstream incom- 
ing plasma approaches the shock transition layer, as seen again in 
the shock front rest frame. This magnetic barrier may reflect back a 
fraction of the incoming protons. In the strongly magnetized case, 
(Ju > 0.1, the coherent gyration of the incoming protons and elec- 
trons gives rise to the emission of large amplitude electromagnetic 
waves through the synchrotron maser instability (Langdon et al. 
1988). As these waves travel back upstream, they lead to a form 
of Wakefield acceleration of the electrons, which leads to electron 
heating at the expense of the incoming protons (Lyubarsky 2006, 
Hoshino 2008). 

The rise of an electrostatic barrier can be described as fol- 
lows. To reach thermaUzation through the shock transition layer, 
the electrons need to absorb part of the wave energy, possibly 
through a kind of anomalous Joules heating. Whatever the mech- 
anism, the length required for the electron heating is several in- 
ertial length 5e = c/wpe, much smaller than the length scale 
of the precursor As electrons remain approximatively in Boltz- 
mann equilibrium, their density increase in the shock transition by 
a compression factor r is accompanied by a potential variation A$ 
such that eA$ ~ Tgln r. If the electrons are heated to equiparti- 
tion with the protons in the downstream plasma, their temperature 
Te ~ Tp ~ -fshi'npc'^ in the downstream plasma. In an electron- 
proton shock, it thus seems unavoidable to have a potential bar- 
rier that reflects part of the incoming protons (see also Gedalin et 
al. 2008); this is an important difference with an electron positron 
plasma. 

If intense waves are excited in the precursor by the reflected 
particles or by emitted electromagnetic waves, their growth may 
also reflect part of the particles, notably electrons. In particular, 
such a mechanism is warranted to ensure the shock transition in 



an unmagnetized pair plasma, as evidenced by various PIC simula- 
tions, see for instance Spitkovsky (2008b). 

The composition of the beam of returning particles is an im- 
portant issue with respect to the development of instabilities in the 
upstream plasma, in particular whether charge neutralization has 
been achieved or not. In that respect, it should be noted that, if elec- 
trons are heated to near equipartition with the incoming protons in 
the shock precursor, the shock itself must behave as if it were a pair 
shock. Then, the electrostatic barrier can be expected to be negligi- 
ble and the returning beam is essentially charge neutralized. 

Whether and to what amount the electrons are heated in the 
shock precursor depends in turn on the instabilities that are gen- 
erated in this precursor. Let us denote by Xe ^ 1 the fraction of 
the incoming (proton) energy carried by the electrons as they reach 
the shock transition layer. Then, Xe = me/ nip if the electrons 
have not been heated in the precursor, and Xe = 1 if equipar- 
tition with the ions has been reached. This parameter Xe is now 
measured by massive PIC simulations (Spitkovsky 2008a, Sironi & 
Spitkovsky 2009, 201 1): generally speaking, it is found that Xe ~ 1 
for strongly magnetized plasmas, Xe < 1 in the intermediate mag- 
netization regime ctu ~ 10"'' — 10~^, and Xe rises up to unity 
again at lower magnetizations. 

As discussed in Lemoine & Pelletier (2010), the returning pro- 
tons and the first generation of accelerated particles form a forward 
beam of energy j^^jUpi? . In the following, we write 7b = 7sh the 
Lorentz factor of the protons composing the beam. If the electrons 
composing the beam have been heated to equipartition with the 
protons before being reinjected towards upstream, either through 
reflection or through shock crossing from downstream toward up- 
stream, they carry a same energy jt\jnpC^ hence their beam plasma 
frequency is also similar. 

In the precursor, the trajectory of a reflected proton is a cycloid 
with an extension upstream that determines the width of the "foot" 
region (i.e. the precursor) 4|sh = "iBhTUpC? / eB-a\a\x (measured in 
the shock front rest frame). As measured in the upstream frame, this 
length scale i!f|u = ^f|sh/7sh; this corresponds to rL|u/7sh with a 
Larmor radius rL|u = 'y%jnpC^ / eB^; the factor 7gij comes from 
the fact that the reflected protons are caught up by the shock front 
after having travelled along the Larmor circle only rL|u/7sh before 
being caught back by the shock front (Gallant Sc Achterberg 1999, 
Achterberg et al. 2001), keeping in mind that the distance between 
the particle and the front is reduced by a factor 1 — ^Sgh — l/(27sh) 
(Milosavljevic & Nakar 2006, Pelletier et al. 2009). 

When the Fermi process develops through cycles of particles 
crossing the front back and forth, the first cycle always occurs and 
the particles participating in this first cycle have an energy compa- 
rable to that of the reflected protons. One can thus assimilate the 
reflected protons to those first cycle particles, and the number of 
particles involved in the Fermi process will only slightly increase 
with further cycles. We will note the ratio of the incoming energy 
density converted into the pressure of these supra-thermal particles: 

?b = ^ , (2) 

and Uh denotes the proper density of the beam particles (i.e., mea- 
sured in the shock front rest frame). This parameter ^b approxi- 
mately corresponds to the ratio of the supra-thermal particle den- 
sity over the upstream particle density measured in this same shock 
rest frame. 
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2.2 Instabilities at a relativistic collisionless shock 

In the case of a current neutralized charged beam, meaning a beam 
carrying positive and negative charges but zero net current along the 
shock normal, the leading instabilities at small magnetization are 
the filamentation mode, the two-stream instabilities, in particular 
the Cerenkov resonance mode with electrostatic modes (or with 
modes of wavenumber parallel to the magnetic field in the case 
of a magnetized plasma) or with Whistler modes (for an electron- 
proton plasma). The filamentation instability takes place at small 
real frequencies TZuj ~ and small parallel wavenumber fcy <^ 
k± relatively to the transverse wavenumbers (parallel and tranverse 
are defined relatively to the shock normal), while the two-stream 
electrostatic instability takes place at resonance Tliu ~ ujp, fcy ~ 
oJp/c in the upstream rest frame, with ujp — {Anny^e^ /nie)^^^ the 
background plasma frequency. 

Such charge driven current neutralized instabilities are of 
course generic in the case of pair shocks. However, even in that 
case, one must expect a transverse current to rise at the tip of 
the precursor, as a consequence of charge splitting in the external 
magnetic field. This current then generates a compensating current 
in the background plasma, which induces a Buneman instability, 
which itself leads to efficient heating of the background electrons 
in the shock precursor This issue will be addressed in Sec.|5] 

In the case of an electron-proton plasma, the beam can be 
current neutralized if, as mentioned above, the incoming upstream 
electrons are preheated to near equipartition with the protons dur- 
ing the crossing of the precursor; then, for all practical matters, 
the shock transition resembles that of a pair shock and the return- 
ing beam carries a vanishing current. Of course, in that case as 
well, one expects the transverse current to emerge in a magnetized 
upstream plasma and give rise to a Buneman type instability, see 
Sec. [5] As mentioned before, such near equipartition has been ob- 
served in PIC simulations in both the high and low magnetization 
limits in oblique collisionless shock waves (Sironi & Spitkovsky 
2011). At high magnetization, electron heating in the precursor is 
due to Wakefield acceleration associated to the ponderomotive force 
of the large amplitude waves emitted by the synchrotron maser in- 
stability (Lyubarsky 2006, Hoshino 2008). At low magnetization, 
electron heating appears to be related to the development of micro- 
instabilities in the shock precursor. 

Such micro-turbulent heating can be understood as follows. 
Consider for simplicity the fully unmagnetized limit and select the 
upstream rest frame; assume further that the shock has reached a 
stationary state. The micro-turbulence is excited on a length scale 
^f|^ that differs from the previous value of the foot length Ij^^, as 
we now neglect any pre-existing magnetic field and assume that 
scattering is dominated by the micro-turbulence (see Milosavlejvic 
& Nakar 2006, Pelletier et al. 2009). The magnitude of this length 
scale is such that the transverse momenta of returning particles dif- 
fuse by an amount (Ap^.x) ~ Pb/7sh during their travel time 
27sh^f|u/c in this precursor. The prefactor corresponds, as 
before, to the difference between the return timescale of the accel- 
erated particles and the precursor light crossing time. For return- 
ing/accelerated particles, the momentum pb ~ 7sh"*pC. Diffusion 
of momenta takes place through scatterings on small (plasma) scale 
electromagnetic fields. Such small scale fields equally contribute to 
heating the upstream plasma electrons, provided that the electric 
wave energy content is not much smaller than its magnetic coun- 
terpart. However, the upstream electrons only experience the length 
scale £f|^ during a light crossing time £{^^/c. Consequently, their 
momentum dispersion amounts to ApJ; ~ rrip(? /2 once the elec- 



trons reach the shock front, which corresponds to equipartition with 
the incoming ions. 

To go one step further, one can easily conceive that the pre- 
heating of the electrons to near equipartition in the precursor pro- 
vides the only means to achieve near neutralization of the returning 
particle current. If indeed electrons are not preheated in the pre- 
cursor but simply overturned (by say, the micro-turbulence), their 
Lorentz factor is increased by a factor 7^1^ but their energy remains 
a factor me/ nip below that of the returning proton beam. Con- 
sequently, their penetration length scale is also a factor rrie/mp 
smaller and for all relevant purposes, the beam can be considered 
as essentially composed of protons. 

In short, a net current along the shock normal in the upstream 
plasma may arise if and only if the electrons have not been heated 
to equipartition with the incoming ions by the time they reach the 
shock front. According to the recent PIC simulations of Sironi & 
Spitkovsky (2011), parallel shock waves offer such an example in 
which the electrons reach the shock transition with an energy that 
is significantly less than that of the protons. It is then natural to ex- 
pect current instabilities, such as the Bell (2004) instability - more 
exactly, its relativistic generalization, see Reville et al. (2006) - to 
develop in the upstream plasma (Lemoine & Pelletier 2010). This 
will be discussed in Sec.|5] 

In the forthcoming Sec. [3]and|4l we concentrate on the fila- 
mentation and two stream instabilites and we discuss their develop- 
ment once the finite angular dispersion of the beam of accelerated 
particles has been taken into account, and considering the possi- 
bility that the electrons of the background (upstream) plasma have 
been preheated through the micro-turbulence to relativistic temper- 
atures. Following the approximation scheme discussed in our pre- 
vious paper (Lemoine & Pelletier 2010), we consider the micro- 
instabilities triggered in the ambient plasma by the returning beam 
over a length scale that is limited by a low background magne- 
tization. In that paper, which discussed the instabilities under the 
assumption of a unidirectional beam and a cold ambient plasma, 
we showed that the Weibel-filamentation instability and the oblique 
two stream instability in particular are very weakly modified by the 
magnetic field during their linear growth, notably because the max- 
imum length scale of the precursor is always much smaller than 
the Larmor radius of the beam. Moreover the magnetization of the 
ambient plasma is always assumed much smaller than unity. Thus, 
in this paper, we neglect the contribution of the magnetic field to 
the dispersion tensor. Of course, such a contribution cannot be ne- 
glected when discussing the development of Whistler waves. How- 
ever we do not consider such an instability in the present paper as 
one of our conclusions is that upstream electrons are heated towards 
equipartition, which quenches the development of Whistler waves. 

3 INSTABILITIES WITH A COLD BACKGROUND 

PLASMA AND FINITE BEAM ANGULAR DISPERSION 

3.1 Beam geometry 

Henceforth, we consider a current neutralized beam composed of 
electrons and protons moving with bulk Lorentz factor 7b in the 
upstream rest frame. In this section, we furthermore assume the 
background plasma to be cold. 

We model the beam with the following axisymmetric wa- 
terbag distribution function: 

/b(u) = -^5{ux — u\\)Q{u\ — u^y — u^^) , (3) 



introducing the velocity variables Ui = pi/(mc), where pi de- 
notes the i— component of the beam momentum. The beam thus 
propagates towards -\-x (parallel direction) and suffers from angu- 
lar dispersion in the perpendicular plane, as characterized by the 
ratio tix/u||. 

In the case of particles returning from or through a rela- 
tivistic shock of Lorentz factor 7sh and propagating in the un- 
shocked plasma, the amount of angular dispersion is known to be 
ui_ ~ W||/7sh- This result is dictated by kinematics: given that 
the shock front is always trailing right behind the accelerated par- 
ticle, once the parallel velocity of this latter drops below /3sh, i-e. 
once the peipendicular velocity of this latter exceeds c/7sh, the 
particle is caught back by the shock wave (Gallant & Achterberg 
1999, Achterberg et al. 2001). As the beam (proton) Lorentz factor 
7b ~ 7sh and 7sh ^ 1, this implies u±_ <^ hence uy ~ 7b 
and uj_ ~ 7sh — ^/yb- This hierarchy allows to calculate the sus- 
ceptibility tensor of the beam by neglecting in a systematic way 
and v?^ in front of u^. Of course, one must pay attention not to ne- 
glect Uy, Uz in the poles of the form to — li. - j3c = to — kiUic/'yh- 
The corresponding beam susceptibility tensor Xi'j is detailed in Ap- 
pendix |A] 

In the final expression for y^j, the angular dispersion en- 
ters through the ratio k±u±c/ R\\, with = 7btJ — k^\u\\c [see 
App.IaI. Therefore its impact on the growth of instabilities is de- 
termined by fcx, 7sh and the nature of the instability which de- 
termines 7?||. At Cerenkov resonance, _R|| ~ "fh^k5, with ujk 
the eigenmode pulsation, and 5 a complex number of modulus 
\5\ <Si 1; at leading order, 5 is a cubic root of unity in the limit 
u±_ — >■ (see Lemoine & Pelletier 2010), hence the growth rate 
Tuj = ujiJiS ~ |^?||/7b| (at resonance). Consequently, the ra- 
tio jfcxitxc/ii|| I ~ k±f3±c/Iuj. If k±f3±c/Iuj > 1, the angu- 
lar dispersion terms dominate over the resonance poles, hence one 
expects growth to be inhibited. Formally, the beam susceptibility 
changes structure in this limit, as discussed in App. |A] Physically, 
one recovers the argument of Akhiezer (1975), that the particles 
travel in the perpendicular direction more than one wavelength of 
the unstable mode on an efolding timescale of the instability, hence 
coherence is lost and growth inhibited, as discussed recently by 
Rabinak et al. (2011). Another important instability in the precur- 
sor of relativistic shocks is the filamentation mode, with TZuj ~ 0, 
fc|| ^ k±. Then k±u±c/ R^^ ~ k±/3±c/Iuj as well. 

Assuming, as we do in App. lAl that the transverse component 
of the wavenumber lies along y, the dispersion relation to be solved 
reads 

12 2,2 , b2\/2 2,2,b2\ 

iuj -OJp-k^ + Xxx^ j " '^P - k\\ + Xyy^ J 

- (fc||fcx + x^j/^')' = 0. (4) 

The susceptibility tensor Xij scales as the ratio squared of the beam 
plasma frequency to the background plasma frequency, i.e. 

(^y^^^nii, (5) 

with ^b defined in Eq. l|2j. In the following, we solve numerically 
the dispersion relation for the growth rates of these various insta- 
bility modes and discuss the inhibition of instabilities due to the 
angular dispersion of the beam. In all numerical calculations be- 
low, we fix ^b = O.I- 
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Figure 1. Contour map of the logio of the l.h..s. of the dispersion relation 
Eq. for a cold background plasma, including the beam contribution. Top 
panel: 7b = 10^, the electrostatic eigenmode includes one branch with zero 
imaginary part and one growing mode (as well as one decaying mode, not 
shown here), as indicated by the darker regions. Bottom panel: 7b = 10^, 
the growing mode has disappeared due to the increased angular dispersion 
of the beam. In both plots, it is assumed '^pb/'^p = O.lme/mp, fcy = 
u)p/{P^c) (resonance condition), kj_ = u)p/c. 



3.2 Oblique two stream instability 

The oblique two stream instability corresponds to the Cerenkov 
resonance of the relativistic beam with the Langmuir modes of 
the background plasma. In the case of a magnetized background 
plasma, the electrostatic modes that are excited propagate along 
the background magnetic field. Since the treatment is similar, we 
focus here on the unmagnetized case. 

To evaluate numerically the growth rate, we proceed as fol- 
lows. We impose the resonance condition fc|| = a;p/(/3bc) and 
look for solutions of the full dispersion relation in the half plane 
(TZu), Xlo > 0), including the beam contribution, for each given 
value of k±. As the beam slightly modifies the real part of the 
root, this latter is slightly displaced from its unperturbed value cjp. 
Figure [T] offers an example of this procedure. It shows the loca- 
tions of the roots of the dispersion relation in the {Tiuj,Xuj) plane 
for 7b = 10^ (top panel) and 7b — 10^ (bottom panel), with 
fcx = ojp/c in both cases. The growing mode is clearly seen in 
the top panel, but absent in the bottom panel in which the condition 
fcx/3xc < lu} is violated due to the smaller value of 7b. 
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Figure 2. Growth rate of the obhque two stream instabihty vs kj_c/u!p for 
various values of 71, : from left to light, 71, = 10^, 10"^, 10*, 10^. The 
cut-offs at high frequency are associated to the inhibition of growth due to 
angular dispersion of the beam; ui'^^/Wp = 0.lm,e/mp, fc|| = ujp/ {/3bc) 
(resonance condition). The dotted line represents the analytical solution for 
the growing mode (see Lemoine & Pelletier 2010). 

To leading order in the growth rate reads (see Lemoine & 
Pelletier 2010) 

Therefore, the condition k±l3±c <C Tiu amounts to 

(7) 

One can thus check that, indeed, for 7sh ~ 10 (corresponding to 
7b — 100), the above condition is violated at k±c/ujp = 1, while 
it is satisfied for 7sh = 300 (7b ~ 10''), in nice agreement with the 
numerical evaluations. 

The effect of the beam temperature on the growth of the two 
stream oblique instability has been discussed previously by Bret et 
al. (2005b) and Bret et al. (2010), although in a slightly different 
context. They show that the two stream instability is rather immune 
to angular dispersion contrary to the filamentation mode, as we find 
here, although the maximum growth rate of the oblique two stream 
mode decreases with increasing dispersion; our findings also match 
these conclusions. In particular, as the angular dispersion increases, 
the maximum growth rate is found at smaller values of k± , with a 
reduced growth rate. 

Figure |2] presents the numerical evaluation of the growth rate 
as a function of k± (solid line) for various values of 7b (hence, 
various values of u±). This figure reveals clearly the wavenumber 
cut-offs that correspond to the finite beam angular dispersion. 

The smallest value of 7ah that allows growth of the oblique 
two stream instability corresponds to setting k± — >■ cjp/7b in 
Eq. ^ above, which leads to 

Note that, at values k± < Up/jh, the oblique two stream mode 
has reduced to the standard parallel two-stream instability. Fur- 
thermore, one must keep in mind that the above condition dictates 
whether growth may occur when angular dispersion is taken into 



Figure 3. Same as Fig. |2] but for the filamentation instability. The calcu- 
lation assumes fcy = 0. Only the growth rate for 7b = 10^ is shown, as 
smaller values of 7b did not lead to exact growing solutions. 

account, yet for growth to occur, other conditions must be met. 
Most notably, the growth timescale must be shorter than the precur- 
sor crossing timescale. This latter condition depends on the degree 
of magnetization of the ambient medium, as discussed in detail in 
Lemoine & Pelletier (2010). 

3.3 Filamentation instability 

The filamentation instability appears at small values of k\\ and TZuj, 
as the solutions to the dispersion relation given in Eq.dD.The effect 
of the beam angular dispersion on the growth rate of the filamen- 
tation instability has been discussed by Bret et al. (2005b), Bret et 
al. (2010), and in the context of relativistic shocks by Lyubarsky 
& Eichler (2006) and recently by Rabinak et al. (2011). This in- 
stability turns out to be extremely sensitive to the beam angular 
dispersion, and as soon as 7b ^ W, the growth is inhibited for all 
values of k±. One important difference relative to the oblique two 
stream instability is the smaller growth rate of the filamentation in- 
stability, all things being equal, which implies that the condition 
k±/3±c > Xuj is more easily satisfied at given values of 7b and k± . 
This is notably illustrated by Fig. |3] which shows the growth rate 
of the filamentation mode for various values of 7b. This numerical 
calculation considers the limit fcy — > and finds the roots of the 
dispersion relation for various values of k±, as before. 

In detail, the growth rate of the filamentation instability in a 
cold background plasma is given at leading order by 

/ fc2 ^2 N 1/2 

therefore, Xcj ^ fc_L/3xc implies 

7sh » ?b ( — ) ik±c « ujp) , (10) 

independently of the value of kj_ as long as fc^c <^ Up. At larger 
values of k±, the right hand side of of Eq. dlOt must be multiplied 
by k±c/ujp hence the condition is more stringent. 

For ^b = 0.1, the numerical solution confirms that the fila- 
mentation instability disappears as soon as 7b < 10*, while it ex- 
ists at wavenumbers shorter than ojp/c for 7b — 10'' ■ This agrees 
well with the above discussion. 
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4 INSTABILITIES IN A HOT BACKGROUND PLASMA 
WITH A FINITE BEAM ANGULAR DISPERSION 

As discussed in some length in Sec. |2l one expects the incom- 
ing electrons to be heated to relativistic energies as they approach 
the shock front. Most notably, this preheating has been observed 
in various configurations in the latest PIC simulations of Sironi 
& Spitkovsky (2011). In particular, at low magnetization, the in- 
coming upstream electrons carry about as much energy as the ions 
when they cross the shock transition. If half of the incoming pro- 
ton kinetic energy density is transferred to the electron component 
and both electrons and protons retain a same bulk Lorentz factor, 
then one can show easily that the electrons are heated to a thermal 
Lorentz factor 7e ~ mp/nie in the electron fluid rest frame. 

We thus consider here a hot electron background plasma, in 
the ultra-relativistic limit; we assume that the protons remain cold. 
We also retain the unmagnetized background plasma approxima- 
tion and assume a ultra-relativistic Maxwellian distribution. 



/e(p) 



(11) 



The longitudinal and transverse (with respect to k) permittivities 
then read (e.g. Silin 1960; see also Hakim & Mangeney 1971, Mel- 
rose 1982, Braaten & Segel 1993, Bergman & Eliasson 2001) 
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with7e = 3//i, /i = meC^ / (ksTe) ■ 

In what follows, the relativistic background plasma frequency 



is written Q,r, 



The full dispersion relation reads 
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Consider now the limit k± <^ ^p/c, keeping in mind that 
resonance implies fcy ~ Qp/c. The dispersion relation then boils 
down to that of the longitudinal mode with 



L . b 



. 



(18) 



We now use the form of the dispersion tensor proposed by Braaten 
& Segel (1993) and set uj = /3pfc||c(l + S), with \5\ < 1 as is 
customary. In this limit, the longitudinal response reduces to 



.^.4+^ln 

where ex = (fc — fc||)/fc|| — 
relation is solved provided 



2 

^pb 
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(19) 



<^ 1. Then the dispersion 



7e 



+ ■ 



(20) 



Due to the smallness of the argument of the log in e , the real part 
of is negative, hence the dispersion relation admits growing so- 
lutions. The argument of the log depends on S, hence the following 
solution 



lui ~ 0.3 



^pb kj 



(21) 



4.1 Two stream instability 

In the fully relativistic regime, the longitudinal mode has a refrac- 
tive index smaller than unity for all values of k, hence the Cerenkov 
resonance with the unperturbed eigenmode can never be fully sat- 
isfied. In all rigor, one should derive the dispersion tensor allow- 
ing for corrections of order in nieC^ /ksT, as one is seeking a 
resonance at a refractive index of order l//3b — 1 + l/(27b)^ 
with a background plasma in which the mean Lorentz factor 7e 
is much smaller than 7b. Such a dispersion tensor has been pro- 
posed by Braaten & Segel (1993). It takes the same form as that 
given in Eqs. il2\ . il3l although all kc must be replaced by kv,, 
where Vt is an effective thermal velocity of electrons in the back- 
ground plasma; hence ~ /3eC. Then, one finds that the re- 
fractive index of the longitudinal mode becomes larger than unity 
(i.e. the dispersion relation crosses the light cone) at some value 
kt — -\/3 (ln(47e) — 2)^''^ iVp/c. For k > k,, resonance is then 
possible. However, extrapolating the form of the dispersion ten- 
sor to the light-like region, one finds that the growth rate of the 
Cerenkov resonant mode with the unperturbed eigenmode is expo- 
nentially suppressed. We thus ignore this branch in the following. 

Once the beam contribution to the dispersion relation is taken 
into account, one finds that an approximate resonance takes place 
for UJ ~ flp, fc|| = uj/PhC for small values k± < flp/c. This 
resonance solves exactly the dispersion relation and leads to growth 
of the modes. Neglecting the effect of angular dispersion, this can 
be understood analytically as follows. 



is valid up to a logarithmic correction (also assuming k± 3> 
f2p/7bc). 

In the opposite limit k± 3> fip/c, one finds two branches: 
the Cerenkov resonance with lj ~ kc, which is exponentially sup- 
pressed as mentioned above and the continuation of the above ap- 
proximate resonance, with cj ~ A:||C ~ ilp. This latter, however is 
nothing but a form of filamentation instability, since it corresponds 
to TZuj <C k±c and fcy <C k±. This branch will therefore be dis- 
cussed in the section that follows [see in particular, Eq. ( I28b 1. 

Let us discuss now the effect of angular dispersion, focussing 
on the Cerenkov resonance at fcxc Q,p — TZu) ~ fcyc. Using 
Eq. | |2U above, one finds that angular dispersion can be neglected, 
i.e. k±P±c <C lid provided 



7.7b » lOOCb ^ — 



(22) 



independently of the value of k±, as long as k± ^ Q,p/c. This 
inequality involves both the thermal Lorentz factor 7e of the back- 
ground plasma and the Lorentz factor of the beam 7b — 7^^. Thus, 
growth can occur provided the temperature of the background 
plasma is sufficiently high; the corresponding threshold tempera- 
ture scales as 7^^^ . 

Figure |4] shows the results of a numerical evaluation of the 
growth rate as a function of k± , for two different background tem- 
peratures, with in each case various values of the beam Lorentz 
factor. The agreement with Eq. ( I22t is found to be quite satisfac- 
tory: for 7e — 30, the instability disappears at values 7b < 10'*, 
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Figure 4. Growth rate of the two stream instabihty in a relativistically hot 
background plasma vs k^c/uj-p {u)-p non relativistic plasma frequency of 
the background plasma) for various values of 7(j at mean Lorentz factor 
7e = 30 (top panel) and 7e = 300 (bottom panel) of the background 
plasma electrons. From left to right: 7^, = 10^ (bottom panel only), 7^ = 
10*, 10^. The calculation assumes kuc = f2p//3b- 



while for 7^ = 300, it is exists down to values of 7b ~ 10 . Note 
that the growth rate peaks at a value fcxc ~ cjp/y^ = Sip, as one 
should expect. 



4.2 Filamentation instability 

The filamentation instability in the hot background plasma can be 
recovered in the limit TZuj — >■ 0, fc||C — >■ 0. We thus write u = iw, 
fc|| — and neglect for the purpose of analytical calculations the 
angular dispersion of the beam. This form of instability has been 
discussed recently in Achterberg & Wiersma (2007). The back- 
ground plasma response can be approximated as 



1 + 3 
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Including the beam contribution, the dispersion relation Eq. il4l 
can be rewritten to leading order as 
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(25) 

We have implicitly assumed k±c 2> w in the above equation. 
Equation i25\ can be solved in the three following limits. 

Iffcxc < w^/^Qp^'-^, then 0,1/ {k_LCw) > k^c^ /w'^ > 1 
andQp/k^c^ 3> 1, so that one obtains 



7 l/3o2/3 
fc_LC < tJpJ, Sip' 



(26) 

The inequality written in parentheses corresponds to the assump- 
tion fcxc < w^/^Slp''^ 

Ifojpf Slp/^ < fcxc < Sip, one rather obtains 



Xlo 



^pb ('^pb'fip^' « fc±c « fip) , (27) 



which corresponds to the standard filamentation growth rate in a 
background plasma at high wavenumbers. The width of the band 
here is governed by (Slp/cjpb)"'"'''' — 'y^^^^ {'me/'mp)~^^'^ . 
It should thus not be much larger than unity. 
Finally, if fcxc 2> Sip, one finds 



Xuj ~ VS^^ (flp « fcxc) 



(28) 



Regarding the effect of angular dispersion, one finds that it can 
be safely neglected provided as before, k±P±c <C lou, or 



1/3 1/2 ^.^ ,^-1/3 

7b > Cb 



-1/3 



(29) 



7 l/3o2/3 

assuming Kxc ^ oj^J^ S2p . 

Figure|5]presents the result of a numerical evaluation of the fil- 
amentation growth rate in a relativistically hot background plasma, 
for various values of 7b, taking into account the angular dispersion. 
In the left panel, the mean background Lorentz factor is 7e = 30 
while in the right panel, 7e = 300. As expected the growth rate 
peaks at fcxc ~ Sip. The angular dispersion shuts off the filamen- 
tation instability at values of k± that scale with 7b: indeed, in the 
range k±c 3> Sip, the condition k±l3±c <C Tuj is much more 
restrictive, as the right hand side of Eq. ( 129) is to be multiplied by 
(fcxc/Slp)^ 

However, contrarily to the results for a cold background 
plasma shown in Fig. [3] the growth rate at small wavenumbers is 
not strongly affected by the angular dispersion. In this way, the rel- 
ativistic temperature of the background sustains the filamentation 
instability at small wavenumbers (relatively to Slp/c) and moder- 
ate values of the beam Lorentz factor. 



4.3 Summary 

Before turning to the current driven instabilities, it may prove use- 
ful to summarize briefly the above results. We have discussed how 
the filamentation and the oblique two stream instabilities are af- 
fected by the finite angular dispersion of the beam and the relativis- 
tic temperature of the electrons of the background plasma. In accor- 
dance with previous works (Lyubarsky & Eichler 2006, Rabinak et 
al. 2011), we find that the filamentation instability is strongly sen- 
sitive to the angular dispersion in a cold background plasma. For 
an electron-proton plasma, a Lorentz factor 7sh ^100 shuts off 
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Figure 5. Same as Fig. |4] for the filamentation instability. Tlie calcula- 
tion assumes fcy = 0. In the upper panel, ■je = 30 while in the lower 
panel, 7e = 300. From left to right: 7b = 10^ (lower panel only), 
7b = 10^, 10*, 10^ 



the filamentation instability in a cold background plasma due to the 
finite angular dispersion, the extent of which varies in inverse pro- 
portion to 7sh. Due to its higher growth rate, the oblique two stream 
instability is much less sensitive to the finite angular dispersion and 
growth may take place at small wavenumbers for 7sh as small as 
10 — 20. As we have stressed, this does not take into account the 
limitation that is imposed by the size of the precursor: for the mode 
to actually grow in a shock precursor, its growth timescale must 
also be smaller than the precursor crossing time. This latter con- 
dition is discussed in detail in Lemoine & Pelletier (2010) and, as 
it depends directly on the level of magnetization of the upstream 
plasma, we do not discuss it here for simplicity. 

Once the electrons are heated to relativistic temperatures, the 
picture becomes susbtantially different. Then, the filamentation 
mode becomes rather insensitive to the angular dispersion at small 
wavenumbers, while the oblique two stream instability is strongly 
inhibited. In both cases, the growth rate of the instability is maximal 
at a wavenumber k± ~ fip/c, which varies as the inverse square 
root of the electron temperature. Conversely, a larger electron tem- 
perature implies a larger spatial scale for the mode of maximum 
growth rate. 



5 CURRENT INSTABILITIES 

Section [2^ has discussed the possibility of current drivent instabil- 
ities at relativistic shocks. In particular, we have argued that a net 
parallel current in the beam of returning particles may emerge at 
oblique shock waves (in the upstream frame) if the electrons of the 
upstream plasma are not heated to equipartition with the protons by 
the time they are overtaken by the shock front. We have also indi- 
cated that, independently of charge neutralization of the returning 
beam, a net peipendicular current may rise through charge splitting 
in an external magnetic field. In the following, we discuss these 
possibilities in turn. 



5.1 Parallel current 

If the shock foot does not preheat the incoming electrons efficiently, 
a parallel current with typical (upstream rest frame) intensity 

ill ~ Cblslwuec (30) 

rises in the shock precursor. Even if a substantial fraction of the 
incoming electrons were reflected at the shock front, the above cur- 
rent would rise, simply because these electrons would carry an en- 
ergy ~ me/nip that of the reflected protons (in the absence of 
preheating), hence their penetration length scale in the upstream 
would be very small compared to that of the protons. 

This current gives rise to a compensating return current in the 
upstream plasma, which may destabilize the pre-existing magnetic 
field on large spatial scales, see Reville et al. (2007) for parallel 
relativistic shock waves. 

It is interesting to note that, for realistic values of and 
7sh, the current exceeds nuCC, hence it cannot be simply com- 
pensated by a non-relativistic drift of the upstream electrons rel- 
atively to the ions. In order to understand how compensation takes 
place, it is instructive to go the shock frame, in which the paral- 
lel current jy |sh ~ Cb7sh'^uec is associated to the charge density 
Cb7sti?iue. Note that the parallel current can be maintained in the 
shock front frame only if there is a net flow of particles escaping to- 
ward upstream; this therefore implicitly means a subluminal (paral- 
lel) shock configuration. Although the shock structure is not strictly 
speaking stationary in this case, we assume here that it evolves on 
long timescales. In this frame, the incoming electrons and protons 
enter the precursor with Lorentz factor 7sh and density 7sh'i-u- If 
the electron fluid is suddenly slowed by an obstacle from velocity 
Psh down to velocity Pex, its shock frame density jumps by a fac- 
tor Psh/Pex but the electron current jex ~ fexPexnexSC - which 
is positive as the electrons flow in the negative x direction - re- 
mains conserved. Therefore, the total incoming e+ p cuiTent must 
remain zero. Given that immediately before entering the precursor, 
the incoming e + p current vanishes, while immediately within the 
precursor, the cosmic ray current does not vanish, it is easy to see 
that one cannot balance the current by simply slowing down the 
incoming electrons or protons. 

In order to achieve charge and current neutralization within 
the precursor, it is actually necessary to reflect part of the electrons. 
Let us assume that a fraction Ke is reflected back towards upstream 
at entrance into the precursor and that electrons and protons move 
with respective velocities jSex and Ppx within the precursor The in- 
coming proton current jpx = — 7sh/3shn.uec, the incoming electron 
current jex ~ 7sh/3sh(l — Ke)nuec, the incoming proton charge 
Ppx = Jpxripxe = jshnuepah/fipx and the incoming electron 

charge p^x = —Jexriexil — Ke)e = — 7sh?^u(l - Kc)el3sh/ Pex- 
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Charge and current neutralization within the precursor thus requires 

Ah 



5b + /3sh/Ppx 



(1 



(31) 



In this case, most of the compensating current in the background 
plasma is carried by the incoming protons, while the repelling and 
the slowing down of the electrons ensure charge balance. As the 
incoming protons are further decelerated by energy exchange with 
the incoming electrons through micro-instabilities, the above rela- 
tion between Pex and j3px must remain valid to preserve charge 
neutralization. Of course, itself varies with location, decreasing 
roughly exponentially ahead of the shock front. Then, the above in- 
dicate that electrons are slowed down progressively as they cross 
the precursor; in the limit /3px — > /3sh, Pex -> (1 — 2^b), leading 
to efficient slowing-down of the electrons. Clearly, one should seek 
this effect in current particle-in-cell simulations. 

The PIC simulations of Sironi & Spitkovsky (2011) indicate 
that the generation of turbulence does not lead to efficient preheat- 
ing of the electrons in the usptream. One may likely relate this to 
the magnetic nature of the modes excited by the Bell instability, 
which do not contain significant electric wave energy in the up- 
stream frame. Nonetheless, the precursor length scale increases in 
time in the case of parallel shock waves, because a fraction of the 
accelerated particles can propagate to upstream infinity in such a 
configuration (Lemoine & Pelletier 2010). Hence, the structure of 
the precursor should itself evolve in time: as the length scale of 
the precursor increases, incoming electrons experience heating for 
a longer duration and arrive at the shock front with a larger frac- 
tion of the incoming ion energy; if this fraction eventually becomes 
comparable to unity, the upstream parallel current disappears and 
the Bell instability shuts off. Unfortunately, present PIC simula- 
tions have not been able to explore the evolution of the precursor 
on such long timescales. 

In the case of oblique shock waves, the simulations of Sironi & 
Spitkovsky (201 1) indicate that the preheating of the electrons does 
not take place efficiently at intermediate magnetisations 10^^ < 
(Tu < 10^^. There is no parallel current in this case in the shock 
frame, only a net charge as discussed in Pelletier et al. (2009). A 
priori, one might expect preheating to occur through a Buneman in- 
stability induced by the perpendicular current (see further below). 
However, as will be shown in the following Sec. 15.21 the magne- 
tization is so large that the Buneman mode does not have time to 
grow before it is advected through the shock front, unless the shock 
Lorentz factor takes quite moderate values. Hence, it should not 
bear a significant impact on the structure of the precursor. 



5.2 Current induced Buneman instability 

Let us assume that the accelerated particle beam carries a net cur- 
rent and induces a compensating return current in the background 
plasma. This situation applies equally well to the case of a paral- 
lel (Sec. 15. il l or a perpendicular current (Sec. 15.31 . The relativis- 
tic motion of the background electrons relative to the background 
ions then induces a relativistic version of the Buneman instability 
(Buneman 1959). 

In the upstream rest frame, the relevant dispersion relation, 
including the relativistic motion of the background electrons with 
velocity vo (Lorentz factor 70) but neglecting the beam contribu- 
tion, takes the form: 



1_ ^ 
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(32) 



The instability develops at frequencies |aj| ^ tJo = k • vo, as the 
small denominator lies in the second term, not the third one. Thus 
we may approximate the dispersion relation as 
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where fctr = [k^ — (k ■ vo)^J represents the component of k 
transverse to the current. The above expression then leads to the 
most unstable mode: 
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As measured in the upstream frame, the advection time across 
the foot is (7ah'^ci)^^, with ojci ~ eBu/{mpc) the proton cy- 
clotron frequency in the background magnetic field. Thus the Bune- 
man instability can effectively grow if T{uik) ^ -yshOJci which 
leads to a condition on the ambient magnetization, as in Lemoine 
& Pelletier (2010): 
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neglecting the 7q''^ dependence. This Buneman mode turns out to 
be the fastest instability, although the oblique two stream mode 
does lie far behind as its growth condition is given by almost 
the same inequality, except that the right hand side term is mul- 
tiplied by . Nevertheless, it is well known (at least in the non- 
relativistic regime) that this Buneman instability saturates rapidly 
by heating the electrons to a temperature such that the anisotropy 
due to the electron current is drowned by the broadened distribu- 
tion (i.e. Vg ~ Ve in non-relativistic regime, being the electron 
thermal velocity, and 7e ~ 7e in relativistic regime), in agreement 
with the discussion of Sec.|3] In practice, and if conditions permit 
it (notably, magnetization), this Buneman mode thus serves as an 
efficient source of electron heating. In turn, this helps the filamen- 
tation mode develop at moderate values of 7sh, which could not 
develop in a cold background plasma given the amount of angular 
dispersion of the beam, as discussed above. 

5.3 Perpendicular current in a magnetic field 

We now consider an oblique shock wave. In the shock front, the 
magnetic field lies perpendicular to the shock normal, to a good 
approximation; we assume that Buish = -Bu|shy. As discussed 
in Sec. 12.11 the returning or acccelerated proton undergoes a cy- 
cloidal trajectory in the Lorentz transformed background field and 
accompanying motional electric field, as measured in the shock 
front frame. Even if the beam does not carry a net parallel current 
at the beginning of the foot, the charge splitting in the background 
field leads to a cosmic ray current oriented along z. Its magnitude 
in the upstream rest frame can be straightforwardly estimated as 
jz.h ~ Cb7sh7^uec, given that the apparent density of cosmic rays 
in this rest frame reads ^b7shnu but that their effective perpendic- 
ular velocity ~ c/7sb. 



5.3.1 Magnetic field structure in the precursor 

The vertical current tends to strongly modify the initial magnetic 
field in a diamagnetic way. Indeed, the choice of orientation of 
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the field, directed toward implies that the vertical current is 
directed toward +2, so that the magnetic field increases toward 
+x, from the shock ramp to the external edge of the foot where 
it reaches its external value. A consistent solution thus requires a 
reduced mean field in the foot. This remains fully compatible with 
the shock crossing conditions, which implies an enhancement of 
the magnetic field strength downstream, when one realizes that a 
similar current develops behind the shock over a distance measured 
by a typical Larmor radius. This current develops in an opposite di- 
rection downstream to that upstream and thus leads to a reduction 
of the magnetic field from its far downstream (asymptotic) value to 
the low value close to the ramp. Note that the cycloidal trajectories 
upstream and downstream have different characteristic gyroradii, 
due to the compression of the magnetic field. 

This diamagnetic effect could be very prohibitive if the per- 
pendicular current were not compensated by the background elec- 
trons. Indeed, in the absence of compensation, one would find an 
induced field Bind, such that 
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(37) 



The modification of the magnetic profile can be calculated as 
follows. First we remark that the motional electric field that com- 
pensates PxBy is uniform throughout the shock transition because 
rotB = 0, hence = Ps\iB^\s\i- Then, in order to characterize 
the compensation of the current, we introduce an effective resistiv- 
ity r) of the background plasma electrons. At the tip of the precur- 
sor, where the beam current is the strongest, one cannot exclude that 
a "double layer" type of structure forms, as in the case of a paral- 
lel shock wave. For the discussion that follows, we ignore this and 
describe the current compensation on phenomenological grounds 
by introducing this effective resistivity. In the precursor and in the 
shock ramp, this resistivity may result from turbulent scattering at 
frequency v^, with r) = A-KV^/iOp^. Ohm's law for the background 
plasma then leads to 



7sh 



(38) 



with j^^pi the background plasma compensating current along z. 
Now the field variation is produced by both the plasma current ,pi 
and the diamagnetic cosmic ray current jz,h'- 



dBy 4tt . , 

OX c 



(39) 



Using Ohm's law to relate jz,pi with By, one thus obtains the fol- 
lowing differential equation that governs the spatial profile of the 
magnetic field: 



£, dBy ^ ^ ^ „ ^ 47r4 . 

7\ Px-Dy — Psh-t>u|sh H Jz,h 

7sh OX 7shC 



(40) 



where £r denotes the resistive length, £r = rjc/iTT = 5^Vs/c. 
The boundary conditions are as follows: for x — > — 00 (far down- 
stream), — >■ —1/3, the current vanish and By — >■ SBuish; for 
X — >■ +00 (far upstream), fix — >■ —fish — —1, the currents vanish 
again and By — >■ B^^^h- In the downstream region in which there is 
a cosmic ray current j^^b < 0, By decreases below 3Bu|sh; corre- 
spondingly, in the upstream region in which jz,h > 0, By increases 
towards its asymptotic value B^^c^i,. 

In the situation that we consider, the magnetic field energy 
is weak and the incoming plasma is able to provide a compensat- 
ing current. The typical length scale over which the compensating 
current is established, £r is much smaller than the length scale over 
which the cosmic ray current is induced ifish, hence one derives the 



modification of the magnetic field atx > £r as By = Buish+AB"*" 
with 

4-Kir 



AB+ ~ 

and the compensating current 

c dB 



47r dx 



J2,b = -(1 H -^jjz.b 

7sh OX 



(41) 



(42) 



The discrepancy with respect to neutralization of the diamagnetic 
current is expressed by the derivative term in the above expres- 
sion, and is of order lr/'ysh£f- Thus the diamagnetic current up- 
stream produces a modification of the magnetic field of a maximum 
amount given by 

-1/2 

• (43) 



AB+ MrUbC _i£r -1/2 

^ 5 ~ 47rCb7sh T-^u 

-Du|sh 7sh-Du|sh "e 



The modification thus remains small if 



^2 I Trie \ -2 



10" 



(44) 



The right hand side is of order W~^'y~^^, hence it should be verified 
in most relevant cases. 

The above transverse current leads to a compensating current 
in the background plasma which induces a Buneman instability, as 
discussed in Sec. 15.21 As the beam transverse current is generated 
at the tip of the precursor, where most of the rotation in the back- 
ground field takes place, the Buneman instability effectively takes 
place in a cold background plasma. It may then lead to efficient 
heating of the electrons, as discussed in Sec. 15.21 



6 DISCUSSION AND CONCLUSIONS 

The present work has discussed the electromagnetic micro- 
instabilities triggered by a beam of shock reflected/accelerated par- 
ticles propagating in the unshocked upstream plasma. In particular, 
it has taken into account the finite angular dispersion of the beam 
of returning particles as well as the possible effects of heating of 
the background plasma electrons. Regarding the development of 
the filamentation and two stream instabilities, the salient results are 
summarized in Sec. 14.31 Let us here discuss how these results af- 
fect our understanding of the development of instabilities at ultra- 
relativistic shock waves, as a function of the shock Lorentz factor 

7sh. 

1) At large values of "/sh ^ 300, both the filamentation and 
modified two stream instabilities can develop in the cold or rel- 
ativistically hot background plasma limits. As mentioned earlier, 
this statement only considers the effect of angular dispersion and 
plasma temperature. For the waves to actually grow, one must sat- 
isfy another condition: the growth timescale must be shorter than 
the precursor crossing timescale. This latter condition depends on 
the magnetization of the upstream plasma and is discussed in detail 
in Lemoine & Pelletier (2010). In the rest of this discussion, we do 
not consider this limitation, which amounts to considering a very 
weakly magnetized upstream plasma. 

2) At lower values of 7sh, the filamentation instability is in- 
hibited by the finite angular dispersion of the beam - the extent 
of which is inversely proportional to j^h ~ in the cold background 
plasma limit, but not in the relativistically hot background plasma 
limit. The oblique two stream instability is inhibited in the hot 
background plasma limit, at least up to some threshold tempera- 
ture which scales as 'y'l^, see Eq. l l22t . 
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3) We have uncovered a situation that leads to the develop- 
ment of a Buneman instability at the tip of the precursor, which 
may efficiently preheat the electrons to relativistic temperatures. 
The Buneman instability is usually triggered by a parallel current. 
But, for the generic case of an oblique relativistic shock wave, the 
rotation of the beam in the background magnetic field leads to a per- 
pendicular current of large intensity. This leads to a compensating 
transverse cuiTent in the background plasma, which in turn leads 
to the Buneman instability, possibly in the relativistic regime. As 
is well known, this instability saturates through the heating of the 
electrons such that the thermal energy becomes comparable with 
the drift energy, thereby drowning the anisotropy of the electron 
distribution. The growth rate of the Buneman instability is larger 
than that of the oblique two stream instability or the filamentation 
mode in the cold background zero angular dispersion limit. 

4) The above then suggests that efficient preheating of the up- 
stream electrons may take place through the Buneman instability 
and through the two stream instability. As soon as the electrons 
have been heated to relativistic temperatures, the filamentation in- 
stability becomes the dominant mode. As indicated by Eq. l |29l ), 
this instability can indeed operate down to values ~ 10 in the 
relativistically hot background plasma limit. For smaller values of 
7sh, the ultra-relativistic shock limit that we have assumed through- 
out is no longer applicable. One should expect different physics to 
come into play in the mildly relativistic limit. 

The upstream electrons are heated to larger thermal Lorentz factors 
7e as they come closer to the shock front, hence the above pic- 
ture must pertain up to the shock front. If the electron temperature 
reaches the threshold discussed previously, the two stream insta- 
bility also becomes efficient. For reference, if the electrons reach 
equipartition with the protons, 7e becomes as large as mp/m^. At 
equipartition, the filamentation instability thus becomes similar to 
that occuring in a pair plasma. As for the Whistler wave instability, 
it develops only if there is a sufficient contrast between the electron 
and proton masses. It thus disappears when the electron relativistic 
mass reaches the proton mass and the modes become right Alfven 
waves, which do not have time to grow unless the precursor has 
been substantially extended through diffusion (see the discussion 
in Lemoine & Pelletier 2010). 

5) The typical perpendicular spatial scale at which the filamen- 
tation and two stream instability growth rates reach their maximum 
scales as ^/T^, due to the evolution of the background plasma fre- 
quency in the ultra-relativistic limit. This suggests that the typical 
spatial scale of the inhomogeneities increases from the electron in- 
ertial scale cjio^c to the proton inertial scale c/tjpi as the modes 
come closer to the shock front. Note that the growth rate of the fil- 
amentation instability does not depend on the temperature of the 
background plasma and takes the same value as in the cold plasma 
limit. 

6) The micro-turbulence that is generated by these instabili- 
ties lead to efficient electron heating. This has been discussed in 
Sec. 12.21 where it has been argued, in particular, that in the unmag- 
netized limit, the large scale of the precursor must guarantee heat- 
ing of the electrons to equipartition with the protons if the wave 
electric energy content is comparable to the magnetic content. Let 
us now account for a finite upstream magnetic field, assuming in 
particular that this magnetic field sets the length scale of the pre- 
cursor. The heating process takes place through the diffusion of 
the electron energy in the bath of micro-turbulence: in fluctuating 
electric fields E on coherence scales l^, electrons reach a thermal 
energy ee such that — 2e^-E^/c^f|u after crossing the precur- 
sor of length The turbulence is supposed to reach a level such 



that it contributes to form the shock, i.e. i5^/(47r) = (^Bn^m^c? 
with being a conversion factor not too far below unity. This im- 
plies a thermal energy ee ~ Cb'^'^'u ^''*7jh^''^'^pC^- For fiducial 
values ~ 10^^, ctu ~ 1Q~^, 7sh ~ 300, the pre-factor is of 
order unity. If one imposes furthermore that the magnetization is 
small enough to guarantee the growth of the filamentation mode, 
one finds CTu < Cb/7sh (Lemoine & Pelletier 2010), the electrons 
are heated to Lorentz factors 7e > ^^/^^^ ^''*?Tip/me, close again 
to equipartition. 

Thus the electrons can be heated in the precursor and roughly ther- 
malized with the protons. 

7) The above considerations agree well with recent PIC simu- 
lations, at least where comparison can be made. In particular, elec- 
tron heating to near equipartition has been observed at small mag- 
netization by Sironi & Spitkovsky (2011), although the detailed 
physical process that is responsible for this heating has not been 
identified in these simulations. The micro-instabilities have been 
observed to take place at moderate values of the shock Lorentz fac- 
tor 7sh ~ 20 and small magnetization (Tu < 10^* (assuming in 
these simulations an electron to proton mass ratio of 1/16). Further- 
more, the typical scale of the fluctuations apparently grows from 
the tip of the precursor to the shock front. 

8) Finally, we question the possible nonlinear saturation of 
the instabilities in the precursor. First of all, one can show that 
these instabilities cannot be saturated by beam particle trapping 
because the timescale to cross a coherence scale is much less 
than the growth timescale of the waves Xu)~^; said otherwise, trap- 
ping of the beam particles would require a prohibitive level of 
turbulence. This can be seen best by going to the rest frame of 
the waves in which the electromagnetic fields are static. In this 
rest frame, the particle cross the transverse coherence length Ic 
in a timescale rni|„ — (Pb|w'c/ei5|w)^^^, which takes into ac- 
count the static force exerted by the fields on the particle in this 
rest frame. As we are interested in the transverse dynamics, one 
can transform the relevant quantities to the upstream frame as 
T"ni|u = 7w|uTni|w, Pb|w ^ Ph ^ 'fsh^pC and E^^ ~ 7w|u-E. Con- 
sider now the example of the filamentation instability, the growth 
rate of which is ^ [m^/nipY^'^LOp. Saturation by trapping would 
require > 7^|u^b7sh('c/5p)^, which is obviously prohibitive. 
This negative conclusion obviously holds equally well for the other 
instabilities that grow faster, in particular the two stream instability, 
the Whistler mode and the Buneman instability. 

Nonlinear effects related to mode coupling thus appear more likely. 
Even when the electrons have turned relativistically hot, the rate of 
nonlinear evolution of the unstable modes ~ as for instance cou- 
pling of oblique two stream modes with transverse waves or acous- 
tic waves - is expected of the order of a;pe£^/(47r7e meC^) ac- 
cording to Zakharov (1972), which leads to saturation once the 
growth rate is balanced by the rate of energy conversion into the 
other stable modes. For the particular example of the two stream in- 
stability, with growth rate ^y^(me/mp)^''^tjp, this leads to ~ 
^i!^ (me/mp)*''^7e, which is therefore not far below ^b if the elec- 
trons have reached equipartition with the ions. As a note of caution, 
one should point out that the above estimates of saturation ignore 
a possible bulk Lorentz factor of the upstream electrons once they 
have been heated in the precursor. Such estimates should neverthe- 
less remain correct in the limit of moderate 7sh. 

The main conclusion of this study of the dispersion and of the 
thermal effects on the growth of micro-instabilities in the foot of 
a relativistic shock, is that electron pre-heating can occur, and that 
this leads to the attenuation of micro-instabilities except the fila- 
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mentation instability, which keeps the same growth rate, the wave- 
length of the instability peak migrating from the electron to the 
proton inertial scale as the electrons are heated to equipartition. 
The plasma that reaches the shock ramp is roughly thermalized 
and behaves like an electron-positron plasma. Hence the following 
question arises: what makes the reflection of a part of the incoming 
particles? Under the assumed conditions of a weak magnetization, 
such that it allows the growth of micro-instabilities, especially the 
filamentation instability, the role of the ponderomotive force ex- 
erted by the growing waves is probably more important than the 
electrostatic barrier This has to be further investigated with dedi- 
cated PIC simulations. 
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APPENDIX A: BEAM SUSCEPTIBILITY TENSOR 

The susceptibility tensor is written as (Melrose 1986): 

Ui d 
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with Ui = pi/(mc), 7 ~ (1 + u^)^'"^ in the above equation and 
rib the beam density (in the upstream plasma rest frame). We recall 
the axisymmetric waterbag distribution function adopted here: 



./b(u) = 



1 



■5{Ux 



(A2) 



The integral in Eq. JAIt can be carried out analytically under the 
approximation discussed in Sec.|3] we neglect in a systematic way 
Uy and Mz in front of but we do not to neglect Uy , Uz in the poles 
of the form uj — k- f3c — lj ~ kiUic/^\^. For simplicity, we rotate 
the perpendicular axes in such a way as to align the perpendicular 
component of the wave vector along y, i.e. kx = k± y. One then 



with 



k±u±c 

7btJ — c , 

U\\UJ — fc||7b c . 



(A 10) 



The beam plasma frequency is defined as usual: cjpb 
[47mbe^/ (7bm)l . The beam velocity /3b — u^^/jh- 
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The above are defined for k±u±c < \R\\\- In the oppo- 
site limit, i.e. when the effects of angular dispersion become sub- 
stantial, the integral in Eg. lAl I contain poles. Using the Plemelj- 
Sohotsky formula to evaluate the integrals over these resonance 
poles, one finds that the above Eqs. ( IA3IA4IA6IA7llA8IA9b for the 
beam susceptibility tensor are continued to the region k±u±c > 
by the substitution (l-x^)^^^ i [x^ - l)^^^ , with 
X = fcxuxc/i?|| in the expressions for the Xij components. 
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